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摘 E: 绿 水 是 维持 黄土 高 原 陆 生 生态 系统 用 水 的 重要 水 源 , 它 与 植被 的 生长 密切 相关 ,分 析 绿 水 


的 变化 特征 并 探讨 影响 其 变化 的 因素 有 助 于 解 » 


干旱 半 干 旱地 区 水 资源 短缺 以 及 生态 环境 安全 
问题 。 以 黄土 高 震 沟 家 区 典型 流域 一 一 砚 瓦 川 流 域 为 研究 对 象 , 采 用 传 抱 开 模 型 评价 了 该 流域 
1981 一 2016 年 的 绿 水 资源 量 ,利用 Mann-Kendall 统 计 检 验 和 滑动 户 检 验 法 对 绿 水 量 进 行 了 趋势 分 
析 和 突变 检验 ,并 用 弹性 系数 法 对 影响 绿 水 量 的 气候 因素 和 下 热 
基于 水 热 耦合 平衡 方程 对 流域 绿 水 变化 进行 了 归 因 分 析 。 结 果 表 明 :(1) 流域 多 年 平均 绿 水 资源 
量 为 503.7 mm, 呈 小 幅 上 升 趋 势 ,2003 年 发 生 突变 ,突变 之 后 绿 水 量 呈 波动 增长 趋势 ;(2) 在 气候 


因素 进行 了 敏感 性 分 析 , 最 后 


影响 因素 中 , 绿 水 对 降水 最 为 敏感 ,潜在 蒸 散 发 次 之 ;(3) 在 气候 变化 和 土地 利用 综合 影响 时 , 绿 水 


对 降水 最 为 敏感 ,其 次 是 下 把 


参数 ,而 对 潜在 燕 散 发 的 敏感 度 最 小 ;(4) 气候 变化 和 土地 利用 变 


化 对 绿 水 变化 的 贡献 率 分 别 为 :74.42% 和 25.58% ,表明 气候 变化 是 导致 绿 水 变化 的 主要 因素 ,其 中 
降水 和 洪 在 蒸 散发 变化 对 绿 水 变化 的 贡献 率 分 别 为 73.63% 和 -1.219%5。 
X 键 词 : 绿 水 ; 气候 变化 ; 土地 利用 变化 ; SHAH, 归 因 分 析 


文章 编号 : 


绿 水 是 指 直接 来 源 于 降水 并 最 终 被 用 于 蒸 散 
消耗 的 水 ”, 是 维持 黄土 高 原 陆地 生态 系统 景观 协调 
和 平衡 的 重要 水 源 , 其 在 支撑 雨 养 农业 和 保障 陆地 生 
态 系统 稳定 发 展 方面 均 具有 不 可 忽视 的 作用 ”。 绿 
水 的 存在 形态 包括 “ 绿 水 流 ”" 和 “ 绿 水 储存 , 绿 水 
流 指 实际 藻 散 发 ,包括 截留 蒸腾、 土壤 和 水 体 表面 
蒸发 ,其 可 以 反映 植被 的 生理 过 程 和 水 文 循 环 过 
程 ,而 绿 水 储存 则 是 存储 在 土壤 中 的 水 分 ,是 支持 
雨 养 农 业 的 直接 水 源 ”。 绿 水 是 水 循环 过 程 的 重要 
组 成 部 分 ,特别 是 在 干旱 半 干 旱地 区 ,其 既 决 定 于 
气候 条 件 , 又 受制 于 土壤 类 型 .植被 种 类 和 土地 利 
用 方式 等 。 因 此 ,在 全 球 气候 剧烈 变化 “以 及 流域 
实施 植被 措施 ”和 工程 措施 ”的 影响 下 ,掌握 绿 水 
的 年 际 变 化 规律 ,探讨 气候 因素 和 土地 利用 因素 的 
驱动 作用 ,可 以 为 干旱 半 干 旱地 区 的 绿 水 资源 高 效 
利用 、 植 被 恢复 等 提供 理论 依据 。 


收 稿 日 期 : 2020-04-26; ”修订 日 期 : 2020-11-13 


关于 绿 水 的 时 空 变化 特征 ,前 人 已 开展 了 一 系 
列 研 究 。 赵 安 周 等 中 基于 SWAT 模 型 评价 了 典型 年 
份 渭河 流域 蓝 水 、 绿 水 的 时 空 差异 ,认为 在 典型 年 
份 绿 水 流 的 变化 相对 于 蓝 水 而 言 较 稳 定 。 吏 婷 婷 
等 运用 SWAT 模 型 估算 了 卢 氏 流域 蓝 水 、 绿 水 资 
源 量 ,分 析 了 蓝 水 、 绿 水 资源 量 的 时 空 分 布 规律 。 
昌 乐 婷 等 “利用 SWAT 模 型 并 结合 SUFI2 算 法 评 
估 了 1977 一 2015 年 及 典型 年 份 细 河 流域 蓝 水 绿 水 
的 时 空 分 布 状况 ,发 现 该 流域 水 资源 构成 以 绿 水 资 
WHE, Gerten 等 ' fi] Hj Lund-Potsdam-Jena (LPJ ) 
模型 评估 了 全 球 尺度 下 土地 利用 变化 对 绿 水 流 的 
影响 ,结果 表明 土地 利用 的 转换 大 大 减少 了 各 个 
地 区 的 绿 水 资源 量 。 综 上 可 知 , 现 有 人 研究 针对 绿 
水 资源 量 进行 评价 时 ,主要 集中 于 对 大 尺度 流域 、 
区 域内 绿 水 资源 量 进行 分 析 与 评价 ,而 针对 干旱 、 
半 干 旱地 区 黄土 高 卉 沟 密 区 典型 小 流域 的 绿 水 变 
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化 还 没有 系统 的 研究 。 此 外 ,在 分 析 绿 水 影响 因素 
时 ,已 有 研究 主要 是 分 析 单 因素 变化 条 件 下 绿 水 
的 啊 应 规律 ,而 定量 分 离 气候 变化 和 土地 利用 变 
化 对 绿 水 资源 量 及 其 时 空 分 布 规律 影响 的 研究 比 
Bk, 

FEF JO, AS SC) HA SC e M AI SRCPCRI ECT f 
域 为 研究 区 域 ,在 对 流域 年 尺度 降水 TEAR UR 
以 及 绿 水 资源 量 进行 突变 检验 的 基础 上 ,定量 突变 
前 后 气候 变化 和 土地 利用 变化 对 流域 绿 水 变化 的 
贡献 率 ,并 探讨 绿 水 资源 量变 化 的 主 控 因素 。 人 研究 
结果 可 对 黄土 高 震 沟 完 区 绿 水 资源 的 可 持续 利用 
和 生态 恢复 等 提供 科学 指导 。 


1 研究 区 概况 


砚 瓦 川 流域 系 泾 河水 系 ,位 于 甘肃 省 庆 阳 市 西峰 
区 及 宁县 境内 (107°37’~107°55’E,35°31’~35°44'N; 
图 1) ,流域 总 面积 为 385.63 km?, 属 于 半 湿 润 到 半 干 
早上 暧 温带 到 中 温带 过 渡 地 带 , 属 于 大 陆 性 季风 气 
候 ,流域 多 年 平均 降雨 量 为 326.8 mm。 流 域 主要 地 
PRTC REI rg FE SB AS 3 PRIS AY , ,是 黄河 水 
利 委员 会 西峰 水 保 站 于 1951 EEZ AY Be e ATARI 
典型 原型 观测 小 流域 ,目前 已 积累 了 大 量 的 降雨 、 
径流 等 实测 数据 。 自 20 世 纪 70 年 代 以 来 , 砚 瓦 川 
流域 治理 经 历 了 4 个 主要 阶段 , 即 1981 一 1993 年 、 
1994—1999 年 .2000 一 2005 年 以 及 2005 一 2010 年 ， 
治理 工作 包括 修建 水 库 和 洲 地 坝 、 人 工 千林、 人工 
种 草 以 及 水 土 保持 生态 工程 等 。 
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图 1 砚 瓦 川 流 域 地理 位 置 及 测 站 分 布 图 
Fig. 1 Geographical location of Yanwachuan Basin 
and the distributions of meteorological, 


precipitation and hydrological stations 


2 数据 来 源 与 研究 方法 


2.1 数据 来 源 

本 文 的 研究 时 段 是 1981 一 2016 年 ,降雨 和 径流 
资料 均 为 实测 资料 。 考 虑 到 资料 系列 的 连续 性 以 
及 雨量 站 的 分 布 情况 ,选用 了 系列 连续 且 均 匀 分 布 
在 流域 内 的 10 个 雨量 站 (图 1) ,通过 泰 森 多 边 形 法 
求 得 流域 的 面 雨 量 。 径 流 资料 来 源 于 流域 出 口 控 
制 站 砚 瓦 川 水 文 站 (图 1)。 气 象 数据 来 源 于 中 国 气 
象 数 据 网 (http://data.cma.cn/) ,包括 平均 气温 、 相 对 
湿度 太阳 辐射 .日照 时 数 和 风速 等 。 

根据 流域 治理 的 状况 ,选用 了 1989 .1995 .2000 
年 和 2010 年 砚 瓦 川 流 域 1 : 100000 的 土地 利用 数 
据 ,结合 流域 DEM 范围 ,利用 ArcGIS 切割 ,得 到 砚 
瓦 川 流域 4 期 典型 年 份 的 土地 利用 数据 ,以 此 为 基 
础 分 析 土 地 利用 变化 对 绿 水 的 影响 。 
2.2 研究 方法 
2.2.1 绿 水 资源 量 计算 模型 ”本 人 研究 采用 傅 抱 瑞 模 
型 "来 计算 砚 瓦 川 流 域 绿 水 资源 量 。 傅 抱 开 模型 " 
是 我 国 气候 学 家 傅 抱 开 教 授 在 Budyko 假 设 " 的 基 
础 上 通过 量 纲 分析 和 数学 推导 出 的 实际 蒸 散 发 计 
算 公式 ,与 其 他 蒸 散 发 计算 公式 相 比 ,该 模型 综合 
考虑 了 下 垫 面条 件 对 绿 水 资源 的 影响 ,其 已 被 证 实 ， 
许多 流域 应 用 效果 良好 "”。 模 型 具体 形式 如 下 : 


Vo 
CI à 
P P P 

式 中 :P 为 降水 量 (mm) ; ET, 为 绿 水 资源 量 (mm); 
ET 是 潜在 蒸 散 量 (mm ) ,本 研究 使 用 Priestley-Taylor 
公式 ”对 其 进行 计算 ;w 为 表征 下 垫 面条 件 的 参数 ， 
其 取 值 范 围 为 (1,c) ”。 
2.2.2 模型 验证 

(1) 水 量 平衡 法 

对 于 一 个 闭合 流域 ,年 尺度 的 绿 水 水 循环 过 程 
可 表示 为 : 


E+AW=P-R (2) 
式 中 :天 +A 和 为 流域 绿 水 资源 量 (mm); P 为 降水 
ft (mm); R HLE (mm), 

本 研究 将 整个 研究 时 段 分 为 2 阶段 ,1981 一 2007 
年 的 数据 用 来 校准 伟 抱 开 模 型 中 的 参数 ,2008 一 
2016 年 的 数据 用 来 验证 参数 ,以 水 量 平衡 法 的 计算 结 
果 为 标准 ,对 传 抱 瑛 模型 中 的 计算 精度 进行 检验 "”。 
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(2) 评价 指标 
通过 试 错 法 ” 率 定 傅 抱 瑛 模型 中 的 参数 w, 用 
与 水 量 平衡 法 计算 结果 的 逐年 相对 误差 ( R. ) 和 均 
方 根 误差 (RMSE ) 检 验 计算 精度 ,并 通过 调整 w ,使 
R, 和 RMSE 最 小 。 计 算 公 式 如 下 : 
R.=(E.,-E,,,) 人 (3) 


RMSE = 125. -Eg (4) 


IP: n 为 研究 时 段 ; Ei, 为 第 ;年 模型 模拟 的 绿 水 资源 
量 (mm); 无。 为 水 量 平衡 法 计算 的 绿 水 资源 量 (mmy)。 
2.2.3 突变 检验 Mann-Kendall(M-K) 突 变 检 验 是 
常用 的 非 参 数 检 验方 法 ,和 常 被 用 于 水 文 气象 长 系 
列 数据 中 各 种 变量 的 突变 检验 ™ ,该 法 可 以 大 致 的 
确定 突变 点 发 生 的 年 份 ,本 文 对 绿 水 资源 量 的 突变 
点 检验 用 M-K 法 进行 分 析 , 并 通过 滑动 下 识别 与 检 
Ja is XE NETT Br WE ,滑动 下 检验 法 既 可 以 对 变异 点 
进行 识别 ,也 可 以 对 变异 点 进行 检验 。 
2.2.4 土地 利用 特征 及 归 因 分 析 

(1) 土地 利用 变化 面积 比例 

土地 利用 变化 面积 比例 3 可 以 反映 流域 所 有 
类 型 土地 变化 总 面积 的 整体 情况 ,表征 流域 土地 利 
用 变化 的 剧烈 程度 。 计 算 公 式 如 下 : 

D= A/S x 100% (5) 

式 中 :D 为 土地 利用 变化 面积 比例 (%); 4 为 所 有 
类 型 土地 变化 总 面积 (km2); S 为 流域 总 面积 (km?)。 

(2) 土地 利用 变化 强度 指数 

土地 利用 变化 强度 指数 能 够 反映 出 某 种 土 
地 利用 类 型 增加 或 流失 的 程度 。 其 计算 公式 如 下 : 

U,-U, 

= COS x 


1 
I y * 100% (6) 


式 中 :7 为 强度 指数 (%); U, 和 U, 表 示人 研究 时 段 初期 
和 末期 该 土地 利用 的 数量 (km); 7 为 研究 时 段 ; 其 
他 参数 意义 同上 。 

(3) 气候 弹性 和 下 垫 面 弹性 系数 的 确定 

本 研究 采用 基于 Budyko 思想 的 弹性 系数 法 
来 分 析 绿 水 对 气候 要 素 和 下 垫 面 参数 的 敏感 性 ,该 
法 量 纲 为 1, 可 以 评价 量 纲 不 同 的 变量 。 根 据 相 关 
研究 ,本 研究 分 别 定义 了 绿 水 的 降水 PRE ZR 
以 及 下 执 面 参数 的 弹性 系数 ,计算 公式 如 下 : 


der 


x (7) 


Ep= 


P 
W, 


L-1 
HG E (8) 
0 g 


I(ET; +P”) , 


eee P tux 


w 


ET? 4 P^ V, (9) 

式 中 : ep 为 绿 水 对 降水 的 弹性 系数 ,表示 单位 降水 
变化 导致 流域 绿 水 量 的 变化 ; ee 为 绿 水 对 潜在 蒸 
散发 的 弹性 系数 ,表示 单位 潜在 蒸 散发 变化 导致 流 
域 绿 水 量 的 变化 ; se。 为 绿 水 对 下 垫 面 参数 的 弹性 
系数 ,表示 单位 下 垫 面 参数 变化 导致 流域 绿 水 量 的 
变化 ; W, 表示 绿 水 资源 量 (mm )。 

(4) 归 因 分 析 

将 研究 期 分 为 突变 点 前 后 2 个 时 段 ,前 后 2 个 
时 上段 多 年 平均 绿 水 量 的 变化 记 为 AW, , AW, 是 由 气 
候 变 化 (降水 ,潜在 蒸 散发 ) 和 土地 利用 变化 引起 
的 ,因此 ,可 将 AW, 表示 为 : 

AW, = AW,» + AW,r + AW,, (10) 

式 中 : AW,, 为 降水 变化 引起 的 绿 水 量变 化 ; AW, 
为 潜在 蒸 散 发 变化 引起 的 绿 水 量 ; AW, AP A 
参数 变化 引起 的 绿 水 量变 化 。 AW,, AW er 、AW,, 
计算 公式 为 : 


ET? ee T 


W, 
AW,» =€ pp AP (11) 
W, 
AW er, = Em pE AET, (12) 
0 
W 
AW. 76, Aw (13) 


式 中 : AP AET, V Ae 分 别 为 突变 点 前 后 2 个 时 段 多 
年 平均 降水 量 潜在 蒸 散 发 以 及 下 垫 面 参数 的 差 值 。 

降水 潜在 茸 散 发 和 下 垫 面 参数 对 绿 水 变化 的 
贡献 率 计 算 公式 为 : 


AW 
= gP 
"AW x 100% (14) 
E 
"e^ ay x 100% (15) 
AW, 
0,7 AW, x 10096 (16) 


RP: tp «Men, m, 分 别 为 降水 潜在 燕 散 发 以 及 下 
热 面 参数 对 绿 水 的 贡献 率 。 
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3 结果 与 分 析 


3.1 流域 气候 特征 分 析 

表 1 是 砚 瓦 川 流 域 基本 水 文 气候 特征 值 ,由 表 1 
可 知 ,1981 一 1989 年 和 1990 一 1999 年 降水 量 减 少 ， 
而 1999 年 之 后 降水 量 呈 现 增 加 的 趋势 ;以 1999 年 
为 分 界 点 ,潜在 蒸 散发 量 呈 现 先 增 大 后 减少 的 趋 
势 ,在 1990—1999 年 潜在 蒸 散 发 增 大 到 峰值 ,之 后 
持续 减少 ;在 整个 研究 期 ,气温 呈现 出 持续 增长 的 
趋势 ,与 全 球 气 候 变 暖 大 背景 一 致 ;干旱 指数 在 
1990—1999 年 达到 最 大 ,最 大 值 为 3.83 ,在 2000 年 
以 后 ,干旱 指数 降低 。 


表 1 砚 瓦 川 流 域 基本 水 文 气候 特征 值 
Tab.1 Basic hydrological and climatic characteristics 


of the Yanwachuan Basin 


年 份 降雨 /mm ”潜在 蒸 散 发 /nm 气温 /CC 干旱 指数 
1981—1989 541.2 1735.0 8.4 3.21 
1990—1999 474.6 1818.8 9.1 3.83 
2000—2009 540.0 1773.5 10.0 3.28 
2010—2016 564.1 1762.2 10.2 3.12 


UE CPGE TBS EK 


从 年 际 变 化 趋势 上 看 (图 2) ,该 流域 降水 量 和 
潜在 蒸 散发 量 呈 现 微弱 的 增加 趋势 ,而 气温 呈现 显 
著 的 增加 趋势 ,于 1996 年 出 现 明 显 谷 值 , 于 2013 年 
出 现 明显 峰值 。 由 图 3 可 知 , 降 雨量 从 1985 年 开始 
降低 ,在 2000 年 前 后 降 至 最 低 , 之 后 持续 增加 ,在 
2000 年 之 后 主要 以 正 距 平 为 主 ;气温 在 1996 年 之 前 
为 偏 冷 时 期 ,之 后 为 偏 暧 时 期 ,分 界 比 较 明显 ;流域 
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图 2 砚 瓦 川 流 域 降水 ,潜在 蒸 散发 .气温 年 际 变化 规律 


Fig. 2 Interannual variation of annual rainfall, potential 


evapotranspiration and temperature in the Yanwachuan Basin 
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图 3 砚 瓦 川 流 域 降水 ,潜在 蒸 散发 .气温 累积 距 乎 值 


Fig.3 Cumulative interannual anomaly values of rainfall, 


potential evapotranspiration and temperature 


in the Yanwachuan Basin 


TEE ZR ALA LTE 1989 EF 2003 4E PAT UE AEE EE Be A 
显 ,1989 年 前 累积 距 平 值 持续 减少 ,1989 一 2003 年 
累积 距 平 值 呈 现 上 升 趋势 ,而 2003 年 之 后 ,其 累积 
WESEL SEU TF AGES AS 
3.2 土地 利用 变化 特征 

统计 不 同时 期 砚 瓦 川 流域 土地 利用 各 地 类 面 
积 及 其 占 比 (图 4) 发 现 , 砚 瓦 川 流 域 土地 利用 类 型 
主要 以 耕地 .草地 为 主 ,耕地 面积 占 比 最 大 , 均 在 
50.0% 以 上 ,其 次 是 草地 , 占 比 在 35.0% 左 右 ; 而 水 
域 面 积 较 少 , 占 比 仅 为 0.1% ,说 明 流 域 河流 水 系 不 

图 5 和 图 6 展示 了 各 土地 利用 变化 面积 比例 和 
各 土地 利用 类 型 增加 或 流失 的 程度 。2000 一 2010 
年 土地 利用 变化 面积 比例 最 大 ,说 明 该 时 期 流域 土 
地 利用 变化 最 为 剧烈 ,林地 增加 程度 最 大 ,草地 流 
失 程 度 最 大 ,这 与 该 时 期 实施 退耕 还 林 工 程 以 及 黄 
河水 土 保 持 生 态 工程 等 措施 有 关 ;1989 一 1995 年 土 
地 利用 变化 面积 比例 最 小 ,说 明 该 时 期 流域 土地 利 
用 相对 稳定 ,变化 面积 比例 相对 较 少 。 林 地 面积 在 
1989—2010 年 加 速 增长 ,分 别 为 0.13% , 0.62% 和 
8.38% ;草地 和 居民 用 地 在 1989—1995 年 都 呈现 增 
大 趋势 ,在 1995 年 之 后 呈现 减少 趋势 ,2000 一 2010 
年 减少 速度 最 快 ,分 别 达 到 了 -8.01% 和 -4.85% ;水 
域 面积 在 整个 研究 期 (1989 一 2010 年 ) 一 直 在 减少 ， 
变化 速度 为 0.-0.03% 和 -0.03% ; 耕地 面积 在 
1989 一 1995 年 和 1995 一 2000 年 减少 ,在 2000 一 
2010 年 增加 ,变化 速度 分 别 为 -1.329% , - 2.0690 和 
4.51%. 
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注 : 饼 图 表示 各 土地 利 


类 型 的 面积 占 比 。 


图 4 不 同时 期 砚 瓦 川 流 域 土地 利用 分 布 及 面积 占 比 


Fig. 4 Distribution of land use and area proportion of Yanwachuan Basin in different periods 
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图 5 不 同时 期 砚 瓦 川 流 域 土 地 利用 面积 变化 比例 
Fig. 5 Proportion of land use area change in the 


Yanwachuan Basin in different periods 
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图 6 砚 瓦 川 流 域 土地 利用 变化 强度 指数 
Fig. 6 Intensity index of land use change in the 


Yanwachuan Basin in different periods 
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3.3 傅 抱 正 模 型 模拟 流域 绿 水 量 

在 校准 期 ,根据 流域 多 年 平均 降水 量 、 径 流量 
以 及 潜在 蒸 散发 量 数据 对 模型 中 的 参数 w 进 行 优化 
CB 7) ,得 到 wow 最 优 取 值 ,从 图 7 可 知 , 当 w 值 偏离 
2.806 时 ,相对 误差 ( R, ) . 均 方 根 误差 (RMSE ) 的 绝 
对 值 在 变 大 ,因此 ,确定 流域 w 值 最 优 为 2.806。 基 
于 最 优 参数 取 值 的 傅 抱 瑛 模型 对 砚 瓦 川 流 域 绿 水 
资源 量 模拟 精度 如 表 2 所 示 。 由 表 2 可 知 , R AM 
RMSE 的 绝对 值 在 校准 期 都 最 小 ,验证 期 的 模拟 精 
度 有 所 降低 ,但 只 范围 也 在 10% 以 内 。 


RMSE/10? mm 


| 1 1 —0.03 
2.800 2.802 2.804 2.806 2.808 2.810 
wo 值 


2.0933 1 


图 7 相对 误差 ( R, ) 均 方 根 误差 (RMSE) 与 w 值 的 关系 
Fig.7 Variations of R, and RMSE with w 


R2 ， 传 抱 瑛 模型 年 绿 水 量 模拟 结果 
Tab.2 Statistics ofthe annual green water of Fu Baopu model 
相对 误差 (及 ) —— 均 方 根 误差 
变化 范围 /% 均值 /% (RMSE)/mm 
率 定 期 (1981 一 2007 年 ) 2.806 -5.22~4.30 0.02 0.0209 
验证 期 (2008 一 2016 年 ) 2.806 -2.83-1.17  -0.92 1.0142 


TE: o 为 表征 下 垫 面条 件 的 参数 ,其 取 值 范围 为 (1,=)。 


图 8 给 出 了 传 抱 开 模 型 校准 期 和 验证 期 的 砚 瓦 
川 流 域 多 年 平均 绿 水 量 ,由 图 8 可 以 看 出 ,对 比 流域 
水 量 平衡 法 与 侍 抱 瑞 模 型 得 到 的 绿 水 资源 量 , 拟 合 
情况 较 好 ,总 体 来 看 模拟 精度 较 高 。 这 说 明和 参数 w 
经 过 优化 后 ,模型 能 够 更 准确 的 反映 出 流域 绿 水 资 
源 量 的 动态 变化 。 
3.4 流域 绿 水 变化 趋势 及 突变 点 检验 

根据 1981—2016 年 砚 瓦 川 流 域 绿 水 量变 化 曲 
线 ( 图 9) 可 知 ,流域 多 年 平均 绿 水 量 为 503.7 mm , 绿 
水 量 年 际 变化 范围 在 288.1~746.1 mm 之 间 ,年 际 变 
化 显著 。 从 变化 趋势 来 看 ,年 绿 水 量 以 0.669 mm-a” 
的 速度 微弱 增加 。 根 据 1981 一 2016 年 砚 瓦 川 流 域 
绿 水 量 累积 距 平分 析 结果 (图 10) 可 知 , 在 2000 年 之 
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图 8 率 定期 与 验证 期 年 实际 绿 水 量 模拟 结果 


Fig. 8 Simulation results of actual annual green water during 


the calibration and verification period 
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图 9 1981 一 2016 年 砚 瓦 川 流 域 绿 水 量变 化 趋势 
Fig.9 Variation trend of green water in the 
Yanwachuan Basin from 1981 to 2016 


300 - 一 BPA m 逐年 距 平 值 


绿 水 资源 量 (累积 ) 距 平 值 /mm 


0L LU 09 102 101 101 1 1 1 
1981 1985 1989 1993 1997 2001 2005 2009 2013 2017 
年 份 


J- 
NS 


图 10 1981 一 2016 年 砚 瓦 川 流域 年 绿 水 量 距 平 值 
Fig. 10 Aomaly of annual green water in the 


Yanwachuan Basin from 1981 to 2016 


前 绿 水 量 呈 波动 减少 趋势 ,在 2000 年 之 后 绿 水 量 呈 
波动 增加 趋势 。 选 用 M-K 检 验 法 对 年 绿 水 资源 量 
的 突变 点 进行 识别 (图 11), 可 以 看 出 ,UF 未 通过 显 
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著 水 平 检验 。 两 曲线 在 置信 区 间 内 有 多 个 交点 , 即 
绿 水 存在 多 个 潜在 的 突变 点 ,为 了 进一步 确定 突变 
点 出 现 的 年 份 ,采用 下 检验 曲线 对 以 上 可 能 出 现 的 
突变 年 份 进行 验证 ,结果 如 图 12 所 示 。 从 图 中 可 以 
看 出 ,突变 点 发 生 时 间 为 2003 年 ,这 与 M-K 检 验 的 
一 个 突变 检验 年 份 吻合 。 以 上 结果 表明 , 砚 瓦 川 流 
域 绿 水 量 在 2003 年 发 生 突变 ,这 主要 是 由 于 , 自 
2000 年 以 来 ,在 黄河 水 土 保持 生态 工程 项 目 支持 
下 ,流域 内 加 大 了 水 土 保持 治理 的 力度 。 
3.5 绿 水 的 气候 及 下 执 面 弹性 系数 
为 了 进一步 分 析 砚 瓦 川 流 域 绿 水 资源 量 对 气候 
变化 和 下 热 面 变化 的 敏感 程度 ,采用 弹性 系数 法 对 
影响 因子 进行 了 敏感 程度 计算 ,具体 为 :将 整个 时 
期 分 为 突变 点 前 后 2 个 阶段 ,突变 点 前 (1981 一 2002 
年 ) 定 义 为 基准 期 ,突变 点 后 (2003 一 2016 年 ) 定 义 为 
变化 期 ,再 分 别 对 突变 点 前 后 2 个 时 期 计算 绿 水 对 
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图 11 砚 瓦 川 流 域 年 平均 绿 水 量 M-K 检 验 结果 


Fig. 11 Results of M-K test for annual green 


water in Yanwachuan Basin 
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图 12 UBL) it ARP Sr 5 FS: HHA 
Fig. 12 Sliding F test of annual average green 


water in the Yanwachuan Basin 


气候 和 下 垫 面 参数 的 弹性 系数 ,结果 如 表 3 所 示 。 

由 表 3 可知, 在 1981 一 2016 年 , 绿 水 对 气候 要 
素 中 降水 和 洪 在 蒸 散发 的 弹性 系数 分 别 为 0.93 和 
0.07 ,对 下 垫 面 参数 的 弹性 系数 为 0.18。 

从 突变 点 前 后 2 个 时 期 来 看 ,流域 下 热 面 参数 
由 1981 一 2002 年 的 2.704, 增 加 到 2003 一 2016 年 的 
3.004, 主要 原因 是 流域 2003 一 2016 年 植被 增加 了 5 ; 
绿 水 对 气象 要 素 的 弹性 系数 中 ,对 降水 的 弹性 系数 
FH 0.93356 8 0.94 ,表明 由 降水 引起 的 绿 水 变化 程 
度 增强 了 ;潜在 蒸 散发 的 弹性 系数 由 0.08 减少 到 
0.07 ,表明 由 潜在 蔡 散 发 引起 的 绿 水 变化 程度 减弱 
了 ; 绿 水 对 下 热 面 参数 的 弹性 系数 由 突变 前 的 0.25 
减少 到 突变 后 的 0.15 ,表明 流域 内 下 垫 面 参数 变化 
对 绿 水 变化 的 影响 减弱 了 。 整 体 而 言 ,2003 一 2016 
年 , 砚 瓦 川 流 域 绿 水 对 降水 的 敏感 性 增加 ,对 潜在 蒸 
散发 的 敏感 性 不 变 , 对 下 垫 面 参数 的 敏感 性 降低 。 
3.6 绿 水 变化 的 归 因 分 析 

表 4 是 绿 水 变化 的 归 因 分 析 结 果 ,可见 ,从 基准 
期 到 变化 期 ,年 均 绿 水 量 增加 了 52.07 mm, Hir, S 
候 变 化 使 绿 水 量 增 加 了 38.75 mm ,对 绿 水 变化 的 贡 
献 率 为 74.42% ,在 气候 要 素 中 ,降水 对 绿 水 增加 的 
贡献 率 为 75.63% ,潜在 蒸 散发 对 绿 水 量 增加 的 贡献 
率 较 小 ,为 -1.21%; 相 比 于 气候 要 素 ,土地 利用 变化 


RI 砚 瓦 川 流域 绿 水 的 气候 和 下 垫 面 弹性 系数 


Tab.3 Elasticity coefficient of green water to climate 


and parameter o in Yanwachuan Basin 

年 份 P [mm ET,/mm w 值 €,/% em/% e,1% 
1981—2002 509.6 1785.6 2.704 0.93 0.08 0.25 
2003—2016 553.9 1756.4 3.004 0.94 0.07 0.15 
1981—2016 526.8 1774.3 2.806 0.93 0.07 0.18 
突变 点 前 后 变化 ”44.3  — 292 0300 - 一 - 
HE: P 为 降水 量 ; ET, DEERE o 值 表示 表征 下 垫 面条 件 的 
参数 ; e, 为 绿 水 对 降水 的 弹性 系数 ; sw 为 绿 水 对 潜在 蒸 散发 的 弹 
性 系数 ; 6, 为 绿 水 对 下 执 面 参数 的 弹性 系数 。 


表 4 砚 瓦 川 流 域 绿 水 变化 的 归 因 分 析 


Tab.4 Attribution of green water changes 


during the two periods 
对 绿 水 变化 的 贡献 率 /% 
AW,/mm AW,,/mm AW xr /mm AW,,/mm ] 绿 水 变化 的 页 献 率 /% 
: | ET, w 值 
52.07 39.38 —0.63 13.32 75.63 -1.21 25.58 


TE: AW, 为 前 后 2 个 时 段 多 年 平均 绿 水 量 的 变化 ; AW, 为 降水 变化 
引起 的 绿 水 量变 化 ; AW on, 为 潜在 蒸 散发 变化 引起 的 绿 水 量 ; AW, 
为 下 垫 面 参数 变化 引起 的 绿 水 量变 化 。 
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使 绿 水 量 增加 了 13.32 mm, 对 绿 水 量变 化 的 贡献 率 
是 25.58% ,由 此 可 见 ,在 砚 瓦 川 流 域 , 气 候 变化 是 绿 
水 量变 化 的 主要 原因 ,而 绿 水 对 降水 的 敏感 性 要 高 
于 潜在 车 散发 ,这 是 因为 绿 水 量 主要 来 源 于 降雨 
量 , 降 雨 越 大 ,流域 蒸发 蒸腾 水 量 越 充分 , 绿 水 量 也 
就 越 大 。 


4 讨论 


气候 变化 和 土地 利用 变化 是 影响 流域 水 文 过 
程 变化 的 2 个 关键 因素 ,量化 两 者 对 绿 水 资源 量 的 
相对 贡献 ,可 以 为 未 来 气候 变化 及 土地 利用 变化 下 
流域 绿 水 资源 的 应 用 和 管理 提供 科学 依据 。 在 全 
球 变 暖 以 及 土地 利用 变化 的 背景 下 ,本 研究 为 了 解 
决 黄土 高 壕沟 餐 区 缺 水 问题 ,提供 了 新 思路 ,分析 
了 绿 水 资源 的 变化 特征 ,并 定量 分 析 了 气候 因素 和 
土地 利用 因素 对 绿 水 资源 的 影响 及 其 对 绿 水 变化 
的 贡献 程度 ,研究 结果 可 对 黄土 高 蚌 沟 餐 区 粮食 生 
产 和 植被 恢复 等 提供 科学 指导 ,同时 可 以 为 黄土 高 
原 地 区 水 资源 短缺 问题 提供 新 的 解决 思路 。 

1981—2016 年 砚 瓦 川 流 域 绿 水 资源 量 呈 现 微 
弱 的 上 升 趋势 ,气候 变化 使 得 绿 水 量 从 基准 期 到 变 
化 期 增加 了 38.75 mm, 对 绿 水 变化 的 贡献 率 是 
74.42% ,其 中 ,降水 使 得 绿 水 增加 了 39.38 mm ,对 绿 
水 变化 的 贡献 率 是 75.63%; 潜 在 燕 散 发 使 得 绿 水 减 
少 了 0.63 mm, 对 绿 水 变化 的 贡献 率 为 -1.21%。 本 
研究 的 结果 与 汉 畅 等 ”运用 SWAT 模 型 在 涟 水 流域 
得 到 的 研究 结果 一 致 ,这 一 结果 的 产生 主要 是 因为 
降水 增加 等 气候 变化 因素 引起 的 ,由 于 绿 水 量 主要 
来 源 于 降雨 量 , 降 雨 越 大 ,流域 蒸发 蒸腾 水 量 越 充 
分 , 绿 水 量 也 就 越 大 。 同 时 ,土地 利用 的 变化 也 会 
致使 绿 水 资源 发 生变 化 ,本 人 研究 中 ,土地 利用 变化 
使 得 绿 水 量 增加 了 13.32 mm ,对 绿 水 变化 的 贡献 率 
为 25.58%, 这 一 结论 与 赵 安 周 等 "基于 SWAT 模 型 
在 渭河 流域 的 研究 结果 一 致 。 连 魏 等 中 在 伊 洛 河 
上 游 流 域 也 得 到 了 相似 的 结论 ,但 荣 琨 等 “在 晋江 
西溪 流域 的 研究 显示 土地 利用 变化 所 引起 的 绿 水 
资源 变化 呈 减 少 趋势 。 导 致 这 一 结果 的 原因 主要 
是 林地 耕地 和 草地 面积 的 减少 。 退 耕 还 林 还 草 工 
程 改 变 了 耕地 林地 和 草地 在 总 面积 中 的 占 比 使 得 
耕地 向 草地 和 林地 转换 ,减少 了 坡 面 径流 ,增加 了 
土壤 含水 量 ,使 得 绿 水 量 增加 ;在 本 研究 中 由 于 退 
耕 还 林 工 程 的 实施 ,林地 面积 增加 , 而 耕地 面积 在 


2000 年 降低 ,在 2010 年 时 林地 面积 和 耕地 面积 相 比 
于 前 一 时 期 都 增加 了 。 耕 地 面积 增加 会 使 农业 灌 
溉 用 水 增加 , 增 大 了 土壤 的 含水 量 以 及 作物 的 蒸 散 
发 ,更 多 的 蓝 水 转化 为 了 绿 水 ,导致 绿 水 增加 。 

在 绿 水 变化 的 驱动 因素 中 ,气候 变化 为 主要 影 
响 因子 (74.42%) ,其 次 为 土地 利用 变化 (25.58% ) ， 
i5 Agee FA Zhang E89 Ze 3t E ie i 2] CX 
南小 河沟 流域 和 Zhao 等 55 在 渭河 流域 得 到 的 研究 
相似 。 本 研究 所 使 用 的 气象 数据 来 源 于 西峰 气象 
站 ,但 近年 来 由 于 城市 不 透水 面积 的 增加 ,导致 观 
测 数据 受气 候 变 化 和 土地 利用 变化 共同 影响 ,使 用 
观测 的 数据 来 模拟 流域 绿 水 量 并 分 析 气 候 变 化 对 
其 的 贡献 率 可 能 不 太 准 确 。 但 是 本 项 研究 中 ,更 多 
的 重点 放 在 了 比较 气候 变化 和 土地 利用 变化 影响 
之 间 的 比较 上 ,而 不 是 偏重 数量 上 。 


5 结论 


本 文 以 砚 瓦 川 流 域 为 研究 区 域 , 基 于 水 热 耘 合 
平衡 方程 对 1981 一 2016 年 流域 的 绿 水 变化 进行 了 
归 因 分 析 ,得 出 以 下 结论 : 

(1) 1981 一 2016 年 流域 内 降水 量 、 气 温和 潜在 
蒸 散 发 量 呈 现 上 升 趋势 ,并 于 2000 年 左右 发 生变 
化 ;在 2000 一 2010 年 土地 利用 变化 最 为 显著 ,林地 
增加 程度 最 大 ,草地 流失 程度 最 大 。 

(2) 1981—2016 年 来 砚 瓦 川 流 域 多 年 平均 绿 水 
量 为 503.7 mm, 呈 上 升 趋势 ,并 在 2003 年 发 生 突 变 。 

(3) 各 气候 要 素 中 , 绿 水 对 降水 和 潜在 蒸 散 发 
的 弹性 系数 分 别 为 0.93 和 0.07; 对 下 热 面 参数 的 弹 
性 系数 为 0.18。 绿 水 对 降水 最 为 敏感 ,其 次 是 下 垫 
面 参 数 ,而 对 洪 在 蒸 散 发 的 敏感 度 最 小 。 

(4) 在 气候 变化 和 土地 利用 变化 共同 作用 下 ， 
绿 水 资源 量 增加 了 52.07 mm ,气候 变化 使 绿 水 量 增 
加 了 38.75 mm, FL oT RAE 74.42% ;土地 利用 变化 
使 绿 水 量 增加 了 13.32 mm ,其 贡献 率 为 25.58%; 气 
候 变 化 是 绿 水 变化 的 主导 因素 。 
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Changing trend of green water and its driving factors in 
the gully region of the Loess Plateau 


WANG Shaona, SONG Xiaoyu', LI Lanjun, LI Huaiyou’, LI Yaolin 
(1. State Key Laboratory of Eco-hydraulics in Northwest Arid Region of China, Xi’ an University of Technology, 
Xi'an 710048, Shaanxi, China; 2. Xifeng Experiment Station of Soil and Water Conservation, 


Yellow River Conservancy Committee, Qingyang 745000, Gansu, China) 


Abstract: The gully region of the Loess Plateau is the main rain-fed agricultural area in western China, which 
has long been considered a fragile area, and water shortage is a critical factor restricting the development of this 
area. Significant changes in climate characteristics and land-use status of the Loess Plateau have taken place with 
the proceeding of global climate fluctuation and regional implementation of soil and water conservation mea- 
sures, which will have a far-reaching effect on green water. Therefore, an in-depth study on the response of green 
water to climate and land-use change is required, which would have remarkable theoretical and practical signifi- 
cance in providing useful indications for constructing an ecological environment construction and the improve- 
ment of the efficient use of green water resources in the Loess Plateau. Green water is estimated using the Bu- 
dyko-based empirical model proposed by Fu Baopu, which excels at the estimation of green water at the basin 
scale, to better understand the response of green water to climate and land-use change from 1981 to 2016 in Yan- 
wachuan Basin, which is a typical small basin in the Loess Plateau gully region. The trend of green water and its 
influencing factors were investigated using the Mann-Kendall test, F-test, and elastic coefficient method. In this 
study, to improve the applicability of the Fu Baopu model, the parameters of the model were calibrated and vali- 
dated using the water-balance model and historical hydrometeorological data (e.g., monitored precipitation data, 
evaporation data, and runoff data). The simulation results show that (1) After parameter calibration, this empirical 
model presents a high simulation accuracy and great applicability in estimating green water over the Yanwachuan 
Basin. (2) The annual mean green water of the study area is 503.7 mm, showing an insignificant increase at the 
temporal scale, and its mutation points occurred in 2003. (3) The elastic coefficients of green water to precipita- 
tion, potential evapotranspiration, and underlying surface are 0.93, 0.07, and 0.18, respectively. It can be seen that 
green water is the most sensitive to precipitation, followed by the underlying surface, which is the least sensitive 
to potential evapotranspiration. (4) The green water increased by 38.75 mm from the base period to the change pe- 
riod due to climate changes, and the green water increased by 13.32 mm from the base period to the change peri- 
od due to land-use changes. The contribution rates of precipitation variation, potential evapotranspiration varia- 
tion, and the underlying surface change to the increase of green water were 75.63%, —1.21%, and 25.58%, respec- 
tively. It was clear that the influence of climate change on increasing green water is stronger than that of the un- 
derlying surface. 


Key words: green water; climate change; land use change; Fu Baopu model; attribution analysis 


